Optical properties of MgH 2 measured in situ in a novel gas cell for 
ellipsometry/spectrophotometry 
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The dielectric properties of a-Mgfi2 are investigated in the photon energy range between 1 and 
6.5 eV. For this purpose, a novel sample configuration and experimental setup are developed that 
allow both optical transmission and ellipsometric measurements of a transparent thin film in equi- 
librium with hydrogen. We show that a-MgEb is a transparent, colour neutral insulator with a band 
gap of 5.6 ± 0.1 eV. It has an intrinsic transparency of about 80% over the whole visible spectrum. 
The dielectric function found in this work confirms very recent band structure calculations using 
the GW approximation by Alford and Chou, 1 As Pd is used as a cap layer we report also the optical 
properties of PdH^ thin films. 

PACS numbers: 77.22.Ch,77.55.+f,78.20.Bh,78.66.Nk 



-a 
c 

o 
o 



> 
oo 

o 

m 
O 
-i— > 



I 

O 
o 



X 



I. INTRODUCTION 

Among metal-hydrides^ the magnesium hydrogen sys- 
tem has always occupied a special place. Magnesium 
reacts reversibly with hydrogen to form MgH2. It is thus 
considered to be one of the most important candidates for 
the reversible storage of hydrogen^ due to its lightweight, 
low cost and high hydrogen storage capacity (7.6 wt.% of 
hydrogen) . In spite of the large number of publications on 
Mg-MgH2 only little is known about the intrinsic physical 
properties of this system. The scarcity of data for MgH K 
is mainly due to the experimental difficulties encountered 
when trying to hydride Mg. 4 Nowadays, a great effort is 
made to improve the hydrogen ab/desorption kinetics by 
making nanocrystallinc Mg 5 and/or adding e.g. transi- 
tion metalsSi^ by ball milling. 

Recent theoretical calculations 9 reproduce that MgH2 
undergoes various phase transitions 10,11 as a function 
of pressure. All theoretical calculations published so 
fajSiiSiiiii (using either the local density approxima- 
tion (LDA) or the generalized gradient approximation 
(GGA)) predict band gaps between 3.1 and 4.2 eV for 
a-MgH2. This is smaller than the few sporadic exper- 
imental values reported until now. Krasko 1 ^ mentions 
a value of 5.16 eV for the band gap from unpublished 
work by Genossar. He and Pongi£ determined in an in- 
direct way using Penn's formulaii an average band gap 
of 5.8 eV. Yamamoto et alJ& report an optical transmis- 
sion spectrum for a thin film in which the transmission 
vanishes at 6.05 eV (205 nm). Apart from that Ellinger 
et ali^ found an index of refraction of 1.95 and 1.96 for 
the ordinary and extraordinary rays at 589.3 nm. The 
dielectric properties have not been studied at all. This 
triggered our interest to study the optical properties of 
MgH 2 in detail. 

Another strong reason for our interest in MgH2 stems 
from metal-hydride switchable mirrors. In 1996 Hui- 



berts et al*& discovered that Y and La thin films change 
reversibly from shiny, metallic films to transparent, in- 
sulating films upon hydrogenation either by changing 
the surrounding hydrogen gas pressure or the electrolytic 
cell potentialiSLS In 1997 Van der Sluis et al^ dis- 
covered that all rare-earth (RE) metals exhibit such a 
switchable behaviour. However, all these materials have 
a characteristic colour in the fully hydrogenated state 
because their band gap is in the visible (E g < 3 eV). 
They showed that alloying with Mg results in colour 
neutral switchable mirrors. This is very important for 
applications in e.g. 'smart' windows. In 2001 Richard- 
son et alH reported that Mg z Ni (z > 2) also features 
reversible switching behaviour upon hydrogenation. In 
all these cases (RE-Mg and Mg^Ni), the band gap is 
shifted to higher energies with increasing magnesium 
concentrationi22i2liSSi2SiSi All these alloys disproportion- 
ate upon hydrogenationi2£i22i22iSi This disproportiona- 
tion is also known for bulk RE-Mgi 31 i 32 The available 
data suggest that the shift of the band gap to higher ener- 
gies is due to the formation of MgH2 which is expected to 
have a large band gap. At the same time the reflectance 
in the low hydrogen phase (when the sample is unloaded) 
increases due to Mg which has a high reflect ion At 
intermediate concentrations the coexistence of Mg and 
MgH2 seems to play an important role in the realization 
of a highly absorbing, black state that is observed during 
loading and unloading of RE-Mg alloysi 3 ^ 3 ^ It may also 
play a role in the black state observed in Mg 2 NiH a; . 36 
Thus, to understand the role of MgH2 in Mg-containing 
switchable mirrors it is essential to determine the optical 
properties of MgH2 thin films. 

In this paper we study thin films of magnesium hydride 
with spectrophotometry and ellipsometry and determine 
the dielectric function and the optical band gap. For this 
purpose we use a special substrate geometry and a novel 
type of optical gas loading cell. 



2 



II. EXPERIMENTAL 

MgH2 thin films are made in two steps. First Pd 
capped metallic Mg films are deposited under UHV con- 
ditions on an appropriate substrate. The Pd cap layer 
is necessary to protect Mg against oxidation and to pro- 
mote hydrogen dissociation and absorption. The films 
are subsequently loaded with hydrogen under high pres- 
sure up to the composition MgEb. 

The hydrogenation of Mg to MgH 2 is, however, not 
straightforward as was shown by Krozer, Kasemo and 
others! 3 ^ 38 ! 39 ! 40 ! 41 Palladium capped Mg films exhibit 
unusual kinetics due to the formation of a blocking MgH2 
layer at the interface between Pd and Mg. The MgH2 
layer prevents H to diffuse 4 ^ to the metallic Mg that is 
still present underneath. The formation of this blocking 
layer can be circumvented by starting hydrogenation at 
relatively low (1 mbar) H2 pressure at a temperature of 
100°C. Magnesium films with thicknesses up to 150 nm 
can be fully transformed to MgH 2 in this way^ 9 ^ 3 - 



A. Film deposition 

Thin, polycrystalline films of Mg and Pd are de- 
posited at room temperature in an UHV MBE system 
with a background pressure of 10~ 9 mbar, using mate- 
rial of typically 99.9 % purity. The magnesium films 
are evaporated from a Knudsen cell and covered with 
a 10 nm thick Pd cap layer. These palladium films 
are deposited from an e-beam evaporation unit. Typi- 
cally, we deposit films simultaneously on a 10x10 mm 2 
glassy carbon substrate for Rutherford backscattering 
spectrometry (RBS), 10x10 mm 2 quartz substrates for 
X-ray diffraction (XRD), resistivity and/or AFM mea- 
surements, and on a quartz substrate (0 42 mm, Heraeus 
Suprasil 1) for optical measurements. 



B. Film characterization 

RBS is used to determine eventual contamination of 
the films. For this glassy carbon substrates are used in 
order to separate the Mg signal from the background sig- 
nal of the substrate. An oxygen contamination between 
0.03 < [0]/[Mg] < 0.085 has been found. 

The thickness of the film is measured with a DekTac 3 
or a Tencor Alpha step 200 mechanical stylus profilome- 
ter. The surface structure, both before and after hy- 
drogen loading, is investigated with a NanoScope III 
atomic force microscope (AFM), operating in tapping 
mode using silicon cantilevers. The scanned areas are 
typically lxl and 5x5 /im 2 from which the root-mean- 
square (RMS) roughness is determined. The thickness 
and roughness values from these techniques are used as 
input parameters in the modeling of the transmission, 
reflection and ellipsometric data (see Sec. IIIII . 



Some samples are contacted ultrasonically with four 
30 /im aluminium wires to monitor the resistivity with 
the Van der Pauw method 4 ^ during loading with hydro- 
gen. 

X-ray experiments are carried out with Cu Ka radi- 
ation in a Rigaku 'Rotaflex' RU 200 or Bruker D8 Dis- 
cover X-ray diffractometer to monitor the transformation 
of hep Mg to rutile MgBi2 in a 6-28 mode. 

C. Optical techniques 

Optical reflection and transmission measurements at 
room temperature (RT) are carried out in a Perkin Elmer 
Lambda 900 spectrophotometer in the range 0.5 < hut < 
6.7 eV (2500 > A > 185 nm). The specular and to- 
tal transmission is recorded while the spectrophotometer 
is purged with argon or nitrogen in order to reduce ab- 
sorption by O2 in the ultraviolet (UV), and H2O in the 
infrared (IR). The quartz substrates (without film) and 
Pd samples are measured in reflection geometry from the 
top side (i.e. the metallic side) at near normal incidence 
(8°) in an absolute reflection unit (using a so-called VN 
geometry). 

Ellipsometry measurements (at RT) in the energy 
range 1.0 to 6.5 eV (1240 > A > 190 nm) are carried 
out in a rotating analyzer variable- angle spectroscopic 
ellipsometer (VASE, J. A. Woollam Co. Inc.), using the 
WVASE32 software program for data acquisition and 
analysis. This instrument measures the ratio of the com- 
plex Fresnel reflection coefficients R of parallel (p) and 
perpendicular (s) polarized light. 45 This ratio defines the 
ellipsometric angles ^>(oj) and A(w) according to 

^ = tan(#(w)) cxp(zA(w)) (1) 
R s 

Three angles of incidence (60, 65 and 70°) are used 
to obtain adequate sensitivity over the whole spectral 
range. Standard deviation and ellipsometric data (\I/, A) 
are recorded at each data point as an average of at least 
100 revolutions, and up to 4000 revolutions of the ana- 
lyzer for the most critical data. 

D. Semi-cylindrical substrate 

As the Pd cap layer on top of the very transparent 
MgBi2 layer is strongly absorbing, ellipsometry cannot 
be carried out from the Pd-side. Thus, ellipsometry 
measurements need to be performed from the 'backside', 
through the substrate. Flat substrates would need to be 
so thick that reflections from the front and backside of 
the substrate can be well separated. With a 2 mm di- 
ameter of the light beam the substrate must be thicker 
than 3 mm. However, at large angles of incidence the 
intensity loss in the light beam is substantial due to re- 
flections at the air/substrate interface. Furthermore, at 
energies close to the limit of the ellipsometer (6-6.5 eV), 
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FIG. 1: Semi-cylindrical quartz glass substrate for ellipsom- 
etry and transmission measurements. The upper figure is a 
side view of the substrate, the one on the bottom a top view. 
A perspective view is given in Fig. The angle of incidence 
of light can be varied between 55 and 75° . On the flat part of 
the substrate three samples with different Mg thicknesses are 
deposited. They are all covered with a 12 nm thick Pd cap 
layer. 



the intensity of the light beam is diminishing quickly. 
These limitations can be avoided with a semi-cylindrical 
substrate. 

This substrate is designed for normal incidence of light 
at the ambient /substrate interface and oblique incidence 
at the internal substrate/MgH2 interface (see Fig.^l. For 
the normal incident approximation to be valid, the diam- 
eter of the semi-cylindrical substrate must be large com- 
pared to the size of the light beam (2 mm) . For practical 
reasons we choose a semi-cylindrical substrate with a di- 
ameter of 42 mm. The top part of the semi-cylindrical 
substrate is cut away parallel to the base surface to en- 
able transmission measurements. This design allows the 
angle of the incident light onto the sample to vary be- 
tween 55° and 75° from the normal. Both the flat and 
the semi-cylindrical substrates are made of quartz glass 
(Heraeus Suprasil 1). This material is transparent deep 
into the UV beyond the limit of our ellipsometer and 
spectrophotometer . 

On the top of the large flat face we deposit, under ex- 
actly the same condition, three films with different Mg 
thicknesses (see Fig. 0. This allows us to analyze com- 
positionally identical films. This method makes the de- 
termination of the dielectric function from ellipsometric 
data more reliable 




FIG. 2: Exploded view of the optical gas loading cell for in 
situ ellipsometry and transmission measurements in a hydro- 
gen environment. In this sketch the cell is shown with the 
semi-cylindrical substrate (see Fig. as sample window at 
the front side. Alternatively, a flat window such as the one at 
the backside of the cell, can be used as sample substrate. 



E. Optical gas loading cell and high pressure 
loading chamber 

In order to measure the optical properties of MgH 2 
and PdHa; in situ in equilibrium with hydrogen at various 
pressures, we designed an experimental setup consisting 
of three parts: i) a substrate/ window (see Fig. [IJ, ii) an 
optical gas loading cell, and iii) a high pressure loading 
chamber. Components ii) and iii) are described below in 
more detail. The setup also includes a special substrate 
holder with a sliding mask for the deposition of the films, 
and a holder to attach the optical gas loading cell to the 
ellipsometer. 

A semi-cylindrical or a flat substrate with the Pd- 
capped MgH2 films deposited on top functions as win- 
dow in the optical gas loading cell (see Fig. - F° r ellip- 
sometry the sample/window is illuminated through the 
substrate to measure the 'backside' of the film. For in 
situ transmission measurements a window is mounted on 
the opposite side of the cell as well. The sample can be 
exposed to a controlled hydrogen gas atmosphere during 
the measurements via two tubes connected to a vacuum 
pump and hydrogen gas cylinder. The cell is designed for 
vacuum, but works also reliably up to a few bar hydrogen 
pressure. 

As already mentioned in the introduction of Sec.[n] hy- 
drogenation of Mg to MgH2 can be successfully achieved 
at moderate temperatures (100°C) by starting at low hy- 
drogen pressure i^Mii Therefore, to be sure that our thin 
films of Mg are completely transformed to MgH2 we start 
loading with a H2 pressure of 1 mbar and increase it in 
steps (within a few hours) up to 100 bar H2. To do 
this the optical gas loading cell is mounted inside a high 
pressure loading chamber (see Fig. [3J. This chamber is 
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FIG. 3: High pressure chamber to load the Mg samples to 
MgHh at a pressure of 100 bar Eb. The optical gas loading 
cell shown in Fig.^is mounted inside the chamber. The vari- 
ous parts are: (1) connection to the vacuum and gas handling 
system, (2) a valve to open and close the high pressure cham- 
ber, (3) and (4) valves to open and close the optical gas cell, 
(5) is the electrical feed-through to enable in situ measure- 
ments of the resisitivity and temperature. 



made of stainless steel and proof-pressurized to 200 bar. 
The design is such that the pressure in the high pressure 
chamber is everywhere the same, and it is not neces- 
sary to expose the sample to the ambient after loading 
at 100 bar H2. It is possible to release hydrogen and 
close the optical gas loading cell with valve (3) when a 
pressure of 1 bar hydrogen or less is reached. With valve 
(2) closed as well, the system is disconnected from the 
gas tubes (1) (while still filled with H2). Afterwards, 
valve (2) is opened to equalize the pressure outside and 
inside the high pressure chamber. Then, the chamber is 
opened, and the valve on top of the optical gas cell (num- 
ber (4)) is closed. The cell can then be disconnected from 
the lid of the high pressure loading chamber while still 
filled with hydrogen. The high pressure loading chamber 
is also equipped with an electrical feed-through (5) with 
several wires to permit measurements of the resistivity 
of a sample during hydrogenation and of the tempera- 
ture inside the chamber with a RhFelOO sensor. During 
hydrogen loading of a Mg film the total chamber can be 
resistively heated up to 100°C. 



III. RESULTS 
A. Sample characterization 

Since Mg is a metal and MgH2 an insulator, the time 
evolution of hydrogenation can be followed in situ in real 
time by monitoring the change of the resisitivity4^ This 
allows us to optimize the hydrogen pressure (starting at 
low pressures and increasing it stepwise to 100 bar H2) in 
such a way that no impenetrable MgH2 layer is formed at 
the interface between Pd and Mg. For practical reasons 
we have mounted an extra sample in the high pressure 
loading chamber for resistivity measurements. 

The resistivity of this as-deposited 150 nm Mg film 
covered with 15 nm of Pd at RT is 6.5 fiilcm (the lit- 
erature value for bulk Mg at 20°C is 4.4 //ficrr4&). The 
reflection of this as-deposited film is high in the visible 
and near-infrared regions (~80%). Both the low resistiv- 
ity and the high reflection indicate the good quality of 
the film. After loading the resistivity reaches 680 /ificm 
under 100 bar H2 at 100°C. Since MgH2 is an insula- 
tor one would at first sight expect a much higher value. 
The moderate resistivity found experimentally is, how- 
ever, due to the metallic Pd cap layer that shortcuts 
the MgH2 layer. Moreover, at a temperature of 100°C 
Mg and Pd may interdiffusc to form a Mg-Pd alloy42^ 
This intermixing has been suggested for Pd capped Y as 
wellfSl and was conclusively shown with photoelectron 
spectroscopy recently^ 2 RBS showed an intermixing of 
Mg and Pd in our films as well. This can be due to 
either alloying, interface roughening or both. The net re- 
sult is that a relatively Pd-rich Pd-Mg alloy is formed on 
top of MgH 2 that absorbs some hydrogen but does not 
become insulating and this causes the shortcut. 

In the as-deposited metallic state, hep Mg has a pref- 
erential growth direction, and only the (002) reflection is 
present in the X-ray diffractogram (see Fig. 01(a)). Load- 
ing a thin Mg film in 1 bar H2 at 100°C does not trans- 
form Mg completely to MgH2 (see Fig. Bfb)). Loading 
at 100 bar and 100°C, on the other hand, left no traces 
of metallic Mg. Only the peaks corresponding to the 
tetragonal structure of the rutile typeiS of ct-MgH2 are 
observed (see Fig. EJc)). Such a preferred growth direc- 
tion is not observed for MgH2 where weak signals from 
the (110), (101) and (200) peaks can be seen. Rocking 
curves around the (002) Mg peak and the (110) MgH2 
peaks show that our samples are polycrystallinc. 

AFM measurements revealed a significant difference 
between the as-deposited Pd covered Mg film and the 
fully hydrogenated films (see Fig. [SJ. Mg expands by 
32% in volume when transforming from hep Mg to rutile 
MgH2. 53 Since the film is clamped by the substrate it 
cannot expand laterally and all the expansion must take 
place out-of-plane. With AFM we indeed noted an in- 
crease in the RMS roughness from 5 nm to 14 nm. It 
can be seen as well that our top layer of Pd is cracked. 
With a mechanical stylus profilometer we found a corre- 
sponding increase of the thickness of the film from 113 to 
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FIG. 4: X-ray diffraction spectra of a 188 nm Mg/10 nm Pd 
film (a) as-deposited, (b) loaded up to 1 bar H2 and (c) of 
a 150 nm Mg/15 nm Pd film loaded up to 100 bar H 2 . The 
large background is due to the quartz substrate. 

162 nm. This 43% increase is larger than the expected 
32% volume expansion because the mechanical stylus has 
a tip radius of 12.5 /mi, and hence cannot probe the deep 
valleys seen on the AFM image. During ellipsometry we 
look through the substrate at the backside of our films 
and not from the top side as with AFM. Nevertheless, we 
found that surface roughening of our Pd top layer needed 
to be taken into account when modeling the ellipsometric 
data. 



B. Transmission and band gap of MgH2 

The optical transmission of Pd capped MgH2 films is 
measured in situ in 1 bar H2 using the gas loading cell 
(see Fig-EJ) - FigureEJshows the total, specular and diffuse 
transmission of a 150 nm thick MgH2 film capped with 
12 nm Pd, loaded at 100 bar H 2 and 100°C. The total 
transmission is measured with the optical gas loading cell 
placed at the entrance port of the integrating sphere in 
the spectrophotometer. In this experiment a flat 3 mm 
thick quartz glass substrate is used. Since we look at 
our film from the substrate side, the MgH2 layer is situ- 




FIG. 5: AFM images of the surface (5x5 fJ-m 2 ) of a 101 nm 
Mg/12 nm Pd (a) as-deposited (b) after loading with hydro- 
gen at 100 bar. The root mean square roughness increases 
from 5 to 14 nm due to the 32% volume expansion accompa- 
nied by the transition from hep Mg to rutile MgH2. 

ated 3 mm away from the port of the integrating sphere. 
The specular transmission is measured with the sample 
in the sample compartment, using the direct detector to 
monitor the signal. The difference between these two sig- 
nals is the diffuse (scattered) transmission. It is probably 
due to the rough surface (see Sec. IIII Al and Fig.[Sfb)) of 
our loaded samples. This diffuse transmission, T<j, has a 
strong wavelength dependence and is proportional to 

T d ex i ex (M 4 , (2) 

where A is the wavelength of light. 54 It is clear from 
Fig- Elk) that the diffuse transmission decreases strongly 
above the band gap as the film starts to absorb light. 
The optical band gap E g , can be estimated from the in- 
tersection of a fitted (fno) 4 curve to the data and an ex- 
trapolation of the flank of the absorption edge. Using 
this so called 'Rayleigh method' we find E s = 5.61 eV for 
this sample and E s = 5.67 eV for a second sample. 

Another estimate for E g can be obtained from the 
absorption edge of the transmission spectra using the 
Lambert-Beer law, T(w) = Tq exp[— a(u))d], where a is 
the absorption coefficient, d the film thickness and To 
contains the transmission of the Pd cap layer and the 
quartz substrate. In the region of the absorption edge 
Tq can be considered as constant in our films (see Figs.[7| 
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FIG. 6: (a) Total (solid line), specular (dashed line) and dif- 
fuse (dashed dotted line) transmission as a function of photon 
energy for a 150 nm thick MgEb film capped with 12 nm Pd 
and loaded at 100° C in 100 bar of hydrogen, (b) Detail of the 
diffuse transmission. A fit using Eq. |5]and an extrapolation 
of the transmission edge are shown. The intersection of these 
two curves gives an estimate of 5.61 eV for the optical band 
gap. 

and 111)11 . The frequency dependence of a near the band 
edge is related to the optical gap through^S*^ 

a(w « ^ —J"-. (3) 

For direct, allowed (forbidden) transitions v = h (u = |) 
and for indirect, allowed (forbidden) transition v = 2 
(u = 3). In amorphous material it has been found that 
v = 2 gives the best results. Combining these equations 
gives 

inrn = inr -c ( ^:^ , (4) 

rthj 

and the constants InToj C, & n d E g are determined from 



a fit to the spectra near the transmission edge in the 
interference-free region. 

Applied to the total transmission this 'Tauc procedure' 
gives a gap of 5.48 ± 0.05 eV using v — 2. It was also 
possible to get a Tauc fit with v = 3 and v = I. How- 
ever, using v = 3 gives values that are too low compared 
to the 'Rayleigh procedure' and we might be fitting an 
interference fringe instead of the absorption edge. For 
v = | the quality of the fit is not as good and we obtain 
a gap of 5.75 ± 0.05 eV from the total transmission. 

Since the diffuse transmission increases rapidly at 
small wavelengths (i.e. with increasing energy), it is 
clearer where the absorption starts in this spectrum than 
in the total transmission. Therefore, we conclude that 
the band gap of MgH 2 is 5.6 ± 0.1 eV. 



C. Ellipsometry 

1. Modeling strategy 

Extracting the dielectric function of a layer from el- 
lipsometric data on samples like ours, which consists of 
several thin layers on a substrate, is a complex task. The 
complicated inversion of the ellipsometric data to the di- 
electric function is, however, greatly simplified if the op- 
tical properties of each individual layer is measured in 
separate experiments. For this reason we adopted the 
following strategy to determine the dielectric function of 
MgH 2 . 

We start by investigating the dielectric properties of 
the quartz substrate. The second step is to evaluate 
the optical properties of the hydrogenated Pd cap layer. 
This is done by investigating a 12 nm thick film of Pd 
on quartz, hydrogenate it in 1 bar H 2 using the optical 
gas loading cell and measure it in the ellipsometer. In a 
third step we study the properties of the hydrogenated 
interface region between Pd and MgH2 carefully since 
earlier experiments with Mg-Pd thin films showed that 
interdiffusion starts already at 100°C42*£2i This is done 
by investigating a 'Pd-Mg' alloy layer consisting of 10 nm 
Mg on a quartz substrate covered with 10 nm Pd. This 
sample is loaded with hydrogen at 100°C and 100 bar 
and measured in the ellipsometer. The final step is to 
measure the total stack (quartz, MgH2, Pd-Mg) and to 
extract the optical properties of MgH2 using the optical 
properties of all the other layers. 

An optical model is defined for each sample and in the 
fitting procedure the difference between the calculated 
(cal) and the measured (exp) A)-values (see Sec. Ill C|) 
are weighted with the experimental standard deviation 
a and fitted with a Levenberg-Marquardt algorithm to 
minimize the mean squared error (MSE) according toS^ 
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where N is the number of (\& , A)-pairs, M is the num- 
ber of fitting parameters, and the indices A and de- 
note data points at different wavelengths and angles. In 
most cases also normal incidence transmission data, T, 
are used to improve the accuracy of the determination of 
the dielectric function^ Then, a third term is included 
in the summation in Eq. 

In the modeling we take into account experimental er- 
rors in incident angle and angular spread due to the sub- 
strate design, and film thickness non-uniformity. It is 
difficult to model both the thickness and the dielectric 
properties simultaneously in ellipsometry^S^I Thus, we 
allow the layer thicknesses to vary only slightly around 
our measured thickness values during fitting. The out- 
put from the modeling consists of the best-fit value of 
the Lorentz-Drude parameters (see Eq. [HJ and their 90% 
confidence intervals. 



2. Optical constants of the glass substrate 

The optical constants of the quartz substrates (both 
flat and semi-cylindrical) and the quartz window used in 
the optical gas loading cell (see Fig. |2J) are determined 
using optical reflection and transmission measurements 
and ellipsometry. This is straightforward, and our results 
match the tabulated values of the refractive index from 
the manufacturer as well as those of quartz glass cited 
in Ref. |6(] (see Fig. Efa)). The extinction coefficients 
being not tabulated in Ref. [6Cj are assumed to be zero. 
Our results on the extinction coefficient show a slight 
absorption near and above 6.5 eV, but still below 10~ 7 
(see Fig. EJb)). The corresponding dielectric function is 
used in the consecutive modeling of the metal hydride 
layers. 



3. Optical properties of PdH x 

A 12 nm thick Pd film, deposited on quartz is exposed 
to 1 bar hydrogen at RT in the optical gas loading cell 
and investigated in the ellipsometer (see Fig. |S| for the 
experimental and fitted data). The Pd hydride, PdHj,, 
that is formed is a strongly absorbing metal. Its dielectric 
function e(uj) = e± + ie2 can be adequately parameterized 
with a Lorentz-Drude (LD) model: 
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FIG. 7: (a) Refractive index, n and (b) extinction coefficient, 
k (b) of the quartz substrates as determined from ellipsometry 
and reflection and transmission measurements. For compari- 
son literature data for the refractive index of quartz glass in 
Ref. 160 are shown as a dot-dashed line. 



where the constant accounts for excitations far above 
6.5 eV; the N Drude terms describe the free-carrier re- 
sponse with u> Pt i the plasma frequency of the i-th Drude 
term and r, the relaxation time; the M Lorentz terms 
represent the effect of interband transitions with fj the 
intensity of the j-th oscillator, ujj its energy, and Tj its 
broadening. The relation between the dielectric function 
and the refractive index, n and extinction coefficient, k, 
is: ei = n 2 — k 2 and e 2 = 2nk. The LD parameters 
obtained for PdHj, in 1 bar H2 are given in Table ||| 

In addition to the sample for ellipsometry, an identical 
sample on a flat substrate is prepared for measurements 
in the spectrophotometer. The transmission and abso- 
lute reflection are determined for this sample in 40 mbar 
hydrogen (4% H2 in Ar). At higher pressures it is not 



8 



TABLE I: Lorentz-Drude parameters and their 90% confidence intervals of a 12 nm thick PdH^ layer in 1 bar H2 obtained 
from ellipsometric data (see Fig. 0. MSE = 7.4, too = 1.266 ± 0.371. All parameters are in eV. 

Jj) Vji (j) ^ \/77 Tj 

(1) 3.389 ± 0.901 0.1892 ± 0.019 

(2) 8.656 ± 2.02 1.775 ± 0.147 



(1) 3.418 ± 0.124 8.588 ± 4.6 5.195 ± 0.691 

(2) 6.878 ± 0.119 11.32 ± 4.63 7.713 ± 0.519 

(3) 9.669 11.96 0.5601 




FIG. 8: Experimental and fitted ellipsometric data, (a) 
and A (b) for a 12 nm thick Pd film in 1 bar H2 used to 
determine the dielectric function of PdHx (see Fig. 1100 . The 
mean square error corresponding to the fit is 7.4. 



FIG. 9: Experimental and fitted transmission (a) and reflec- 
tion (b) data for a 12 nm thick Pd film in 40 mbar H 2 (4% H 
in Ar) used to determine the dielectric function of PdH^ (see 
Fig. I10II . The mean square error corresponding to the fit is 
1.04. 



possible to use the optical gas loading cell when we deter- 
mine the absolute reflection since we must look directly 
at the Pd sample and not via the substrate. To obtain a 
hydrogenated Pd sample, the whole spectrophotometer is 
purged in Ar containing 4% H2 , corresponding to a par- 
tial pressure of 40 mbar H2. This is the highest H2 con- 
centration we can use in the (open) spectrophotometer. 
In Fig. the experimental and fitted data are displayed, 
the LD parameters obtained for Pdl^ in 40 mbar H2 are 
given in Table ITT1 

The difference in the MSE (see Table|l|and|nj) between 



the 1 bar and 40 mbar measurements is mainly due to 
the different substrate geometry, semi-cylindrical vs. flat 
substrate. We assign the major part of the difference in 
the MSE to the cylindrical incident and exit surface of 
the semi-cylindrical substrate and to the fact that the 
ellipsomety measurements have been performed from the 
'backside' of the sample through the substrate. Note that 
this does not change the dielectric function, it merely 
gives a larger spread in the input data for the analysis. 

In Fig. we show the resulting dielectric function at 
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TABLE II: Lorentz-Drude parameters and their 90% confidence intervals of a 12 nm thick PdHz layer in 40 mbar Eb (4% H2 
in Ar) obtained from reflection and transmission data (see Fig. EJ. MSE = 1.04, too = 1.280 ± 0.241. All parameters are in 
eV. 

Jj) Vji U) ^ \/7± £j 

(1) 5.156 ± 0.955 0.001 ± 0.0241 

(2) 8.382 ± 2.65 2.238 ± 0.323 

(1) 4.131 ± 0.111 9.815 ± 4.17 7.104 ± 0.581 

(2) 7.623 ± 0.249 7.793 ± 4.14 0.927 ± 0.235 



(a) 




derived from the Drudc parameters using 



1 



3 4 5 

Energy (eV) 



FIG. 10: Real (a) and imaginary (b) parts of e(w) for the 
PdHU films as determined from ellipsometry in 1 bar H2 
(see Fig. |SJ and reflection and transmission measurements in 
40 mbar H 2 (4% H 2 in Ar) (see Fig. OH . For comparison, the 
dielectric function found by Rottkay et aim - for PdlrL in 4% 
H2 and the dielectric function of pure Pd according to Ref. 
are included. 



1 bar, and 40 mbar H2 partial hydrogen pressure, and 
compare it to literature data for PdH^ by Rottkay et 
al£^ and literature data for Pd from Ref. l60L 

The plasma frequency, lu p , of PdH x in 1 bar H2 is 
slightly larger than the one in 40 mbar H2. Because 
Ljp oc tJTQ with n c the charge carrier density, there 
are more free charge carriers in PdH,; in 1 bar than in 
40 mbar H2. Values for the (optical) resistivity p can be 



(7) 



with eo the vacuum permittivity, lu p the plasma fre- 
quency and r the electron relaxation time. This opti- 
cal resistivity calculated for the dominant Drude term 
no. 2 shows the same trend as the plasma frequency it- 
self: In 40 mbar H2 PdH^ has a resistivity of 235 /lilcm, 
in 1 bar 178 /ificm. However, bulk Pd has a resistivity 
of 10.53 /ificm and PdH^ has a maximum resisitivity of 
about 20 /iSlcm when x = 0.7 at RTiS Our much larger 
resistivities are probably due to the fact that the 12 nm 
thick Pd film consists of somewhat disconnected islands. 
Hydrogen absorption causes the Pd islands to expand 
which decreases the resistivity between them&Si. Thus, 
we have a sort of percolation effect and the resistivity is 
lower in 1 bar H2 than in 40 mbar H2 contrary to bulk 



4- Optical constants of the double layer Pd/Mg 

To investigate the optical properties of the partially 
interdiffused Pd-Mg top layer, we deposit a layer of 10 nm 
Mg capped with 10 nm Pd on quartz. This Pd-Mg film is 
then exposed to hydrogen at 100°C and 100 bar (together 
with the thick Mg film covered by Pd (see Fig.[Q)). After 
hydrogenation the optical properties are determined in 
the optical gas loading cell at RT. This sample is treated 
as consisting of two layers: a Mg-rich Pd-Mg alloy on the 
substrate covered with a Pd-rich Pd-Mg alloy on top. 
As starting values for the fitting procedure we use the 
PdHa; dielectric function determined above for the Pd- 
rich Pd-Mg top layer, and combine it with voids in a 
Bruggeman effective medium approximation (EMA) to 
take surface roughness into account^ A Lorentz-Drude 
model is used for the second layer, the Mg-rich Pd-Mg 
alloy. Ellipsometric data for three angles of incidence (55, 
60 and 65°) and normal incidence transmisson data are 
then combined in a multiple data type fit. All data are 
measured on the same sample and during fitting the layer 
optical functions and thicknesses are coupled. The final 
iteration results in a fit with a MSE of 8.7. The optical 
properties of this double layer are then used as starting 
values for the top layer of the thicker MgH2 film. 
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5. Dielectric function of Mgli-z 



In the evaluation of the optical properties of MgH2 we 
analyze ellipsometric and transmission data of a 124 nm 
thick MgH2 film capped with 12 nm PdH^ (as mea- 
sured with a mechanical stylus profilometer in the hy- 
drogenated state). In addition to these data, transmis- 
sion data of a compositionally identical film, but with a 
thickness of 162 nm (when hydrogenated) are included 
in the modeling. These three data sets are evaluated in 
three parallel, coupled models simultaneously. The main 
features can be modeled using two Lorentz oscillators at 
the high energy side of the measured spectra, at 6.4 and 
6.9 eV. These oscillators mark the beginning of the con- 
duction band. 

The optical parameters of the capping layer, consisting 
of the Pd-rich Pd-Mg alloy on top of the Mg-rich Pd-Mg 
alloy, are initially fixed to the parameter values obtained 
in Section IlII C 41 The only parameters of the top layers 
which arc allowed to vary are the thicknesses since the 
diffusion of Pd into MgH2 may be larger than the 10 nm 
in the thin Pd/Mg double layer. In the final iteration a 
global fit is used in which all LD parameters are allowed 
to change. The final MSE is 17.17. Table Ell gives the 
LD parameter values from the final iteration for MgH 2 , 
the Pd-rich Pd-Mg cap layer and the Mg-rich Pd-Mg cap 
layer. 

The plasma frequencies obtained for the two top layers 
give us a clue about their composition. Since the plasma 
frequency of the top layer (lj p — 14.35 eV) is much larger 
than the one of the lower cap layer (u> p = 6.672 eV), the 
top layer has a larger charge carrier density and is thus 
more metallic than the lower one. This indicates that the 
top layer is formed by a metallic Pd-Mg alloy. The lower 
layer contains some insulating MgH^ as well. The optical 
resisitivity (see Eq. 0) of the top layer is 38 [iflcm com- 
pared to an electrical resistivity of 63 /iSlcm (at RT) for 
the total stack as measured after loading. Thus, the top 
layer is indeed shunting the resisitivity measurements of 
MgH2. Since PdH^ has an optical resistivity of 178 /j,Slcm 
it is clear that the top layer contains some metallic Mg 
as well which has a much lower resistivity (6.5 /iSlcm). 

The total thickness of the stack obtained from ellip- 
sometry is 130.6 nm after hydrogenation. With the stylus 
profilometer we found a thickness of 136 nm. However, 
as mentioned before the profilometer gives a value that is 
too large. Before hydrogenation the thickness was 95 nm. 
This would mean an increase of 37.5% instead of the the- 
oretical 32% volume expansion. 

In Fig. Illf ahfb') the experimental and fitted values of 
iff and A are given, in (c),(d) the experimental and fitted 
transmission curves are shown. Finally, Fig. I12f ahfb) 
shows the real and imaginary part of e(ui) obtained for 
cv-MgH 2 . 




Energy (eV) Energy (eV) 



FIG. 11: Experimental and fitted data for a 124 nm thick 
MgEb film covered with 12 nm Pd used to determine the di- 
electric function of MgH2 (see Fig. 1121 ; (a) ellipsometric data 
for if 1 , (b) ellipsometric data for A, (c) optical transmission (d) 
transmission of a sample consisting of 150 nm MgH2/12 nm 
Pd. All experimental data are modelled simultaneously in the 
fitting procedure as described in Sec. 1111 CI The mean square 
error of the fit is 17.17. 



IV. DISCUSSION 

In Table llVl both experimental and theoretical values 
for the band gap of MgH2 are given. Hartree-Fock (HF) 
calculations are not included since this method overesti- 
mates the band gap considerably. Very recent and not 
yet published theoretical work by Herzig 66 , AuluckSi and 
Alford and Choui on MgH2 is included. 

The band gap of 5.6 ± 0.1 eV determined in this work 
for a-MgH2 is close to values mentioned sporadically in 
literature. The value 5.16 eV was obtained in an UV- 
absorption study mentioned by Krasko. 1 ^ However, de- 
tails about how this value was obtained have never been 
published. The value found by He and Pongi& for the 
average band gap < E g > is close to ours. This is rather 
surprising since < E g > was obtained in an indirect way 
from X-ray photoelectron spectroscopy (XPS) data us- 
ing Penn's formula^ Yamamoto et alm& measured the 
specular optical transmission of thin layers of MgH2 cov- 
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TABLE III: Lorentz-Drude parameters and their 90% confidence intervals of a 124 nm Mgtb film covered with 12 nm Pd 



obtained from ellipsometric and transmission data (see 
eV. 


Fig 


ITT). 


MSE = 17.17 and £oo 


= 1.595 ± 0.092. 


All parameters are in 


MgH 2 film, d = 115.9 nm 


(i) w P ,i l/n 


(?) 












(1) 
(2) 




6.4 
6.9 


5.516 

7.689 ± 0.519 


0.6454 

0.01223 ± 0.109 


Top cap layer: Pd-rich Pd-Mg alloy (with 53% voids), 


d = 


15.6 


nm 








(i) 










(1) 14.35 ± 0.062 1.055 ± 30.6 


(i) 

(2) 




2.964 ± 0.122 
8.5 


5.208 ± 0.559 
17.408 ± 12.2 


1.702 ± 9.03 
0.01366 ± 4.96 


Lower cap layer: Mg-rich Pd-Mg alloy, d = 4.1 nm 




U) 










(1) 6.672 ± 0.865 0.003629 ± 8.29 


(i) 




3.293 ± 0.0824 


24.61 ± 0.420 


14.09 ± 7.06 




This work 
Alford and Choir 




Energy (eV) 



FIG. 12: Real (a) and imaginary (b) parts of t(uj) for MgH2 
determined from ellipsometry and transmission data (see 
Fig. [TTJ. The real and imaginary parts calculated by Al- 
ford and Chou— using the GW approximation are shown for 
comparison. The transmission of a 100 nm MgHh film in vac- 
uum calculated with the dielectric functions shown in (a) and 
(b) is displayed in the lower panel (c). The theoretical gap is 
approximately 0.5 eV too large. 



TABLE IV: Literature data for the band gap, E g , of MgHb 
both experimental (expt.) and calculated (calc). In paren- 
theses the smallest direct gap for MgHb is given. 



Material Method E g (cV) Ref. 

a-MgEb expt. UV-absorption 5.16 1^ 

XPS 5.8 ,16 

EELS >3.7 a 68 

XPS >3.4 a 69 

Transmission 6.05 L^ 

Ellipsometry/ 5.6 ± 0.1 this work 
Transmission 

calc. LDA 3.06 12 

LDA 3.4 13 

LDA 3.45 (4.34) 66 

LDA 3.3 67 

LDA 3.10 1, 

GGA 3.78 14 

GGA 4.2 2 

sX-LDA 5.71 (6.41) 66 

GWA 5.25 (6.11) J. 

7-MgH 2 calc. GGA 4.3 9, 

/3-MgH 2 calc. APW 0.23 70 

GGA 2.35 71 



"If there are no charging effects, the band gap would be twice the 
indicated value. 
''This is the photon energy where the transmission vanishes. 



ered by Pd and found that the transmission is zero at 
6.05 eV. However, they did not apply Tauc's method (see 
Sec. IIII 5|) to the transmission edge in order to obtain an 
estimate for the band gap. Furthermore, one should keep 
in mind that 6.02 eV is at the detection limit of their Shi- 
madzu spectrophotometer. 
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As can be seen in Table IIVI LDA calculations give 
a band gap that is systematically too low. This is a 
well-known feature of this approximation. Similarly, the 
GGA used in two other papers2ii4 to calculate the den- 
sity of states of MgH2, underestimates the band gap. Our 
band gap is closest to the theoretically calculated gaps 
of Herzig 66 using screcncd-exchange-LDA (sX-LDA) and 
Alford and Choui using the GW approximation (GWA). 
Furthermore, it is interesting to point out that our ex- 
perimentally found gap E g = 5.6 ± 0.1 cV is very close 
to the difference in ionization energy between Mg and H: 
5.952 eV. 

Back in 1955 Ellinger et alJ^ determined the refractive 
index of a-MgH 2 at 589.3 nm (2.107 eV) and found n = 
1.95 and 1.96 for the ordinary and extraordinary rays, 
respectively. We find n = 1.94 and k = 7.6T0~ 3 at the 
same energy which is very close. 

Both AuluckSL using LDA and Alford and Chou 1 using 
LDA and GWA have calculated the band structure and 
dielectric function for a-MgH2. To obtain the dielectric 
function only direct transitions are taken into account. 
The only difference between the LDA and GWA curves 
is the energy position. The dielectric function obtained 
for a-MgH2 with GWA agrees quite well with our mea- 
sured values (see Fig. ^](a),(b)), while the LDA curve, 
as expected, is shifted to too low energies. This indicates 
that a scissors-operation that shifts the conduction band 
rigidly with respect to the valence band works well to 
correct LDA calculations^^ The energy at which both 
ei and €2 exhibit a marked increase is, however, slightly 
different for the experiment and the GWA calculation. 
sX-LDA seems to overestimate the band gap since it gives 
even larger values for the (in) direct gap than GWA. 

Figure 112( c) shows the optical transmission for a 
100 nm thick MgH2 film in vacuum as calculated with 
our experimental dielectric function and the one calcu- 
lated by Alford and Chou. As can be seen MgH2 has 
an intrinsic transparency of about 80% over the entire 
visible spectrum. The difference in energy between the 
absorption edges of the two transmission spectra and the 
dielectric functions is about 0.5 eV. Thus, GWA overes- 
timates the optical gap by 0.5 eV. The direct gap deter- 
mined from the band structure using GWA is 6.11 eV. 
Subtracting 0.5 eV from 6.11 eV gives a value of 5.61 eV 
which is within the error margin of our experimentally 
found gap. More information about the dielectric func- 
tion and band structure of MgH2 by Alford and Chou 
will be published elsewhere. 

Yet unpublished calculations indicate that the opti- 
cal properties of both a- and 7-MgH2 are close to each 
other>iiSi Vajeeston et al~ find a gap of 4.2 and 4.3 eV 
for a- and 7-MgH2, respectively with GGA; Bastide et 
alMi and Bortz et ((/.— have found that the structures of 
the two different phases are closely related and that the 
density is almost the same as well as the H-H distances. 
Therefore, our dielectric function for a-MgH2 is probably 
a very good approximation for 7-MgH2 as well. This is 
important to model the optical properties of switchable 



TABLE V: Literature data for the band gap, E g , of materi- 
als closely related to MgH2 both experimental (expt.) and 
calculated (calc). 



Material 




Method 


E g (eV) 


Ref. 


MgF 2 


expt. 


Reflectance 


12.4 


M 




calc. 


tight-binding, 










pseudopotentials 


12.8 


75 


MgO 


expt. 


Reflectance 


7.77 


76 




calc. 


LDA 


~5 


II 






GWA 


7.8 


77.78 


MgS 


calc. 


LDA 


2.6 


79 






LDA (corrected) 


4.59 


80 


MgSe 


expt. 




5.6 


81 


MgTe 


expt. 




4.7, 3.6 


81 


CaH 2 


expt. 


XPS 


~5" 


£2 






SRPES 


~5 a 


83 




calc. 


LDA 


3.32 


84 


SrH 2 


expt. 


XPS 


~5 a 


82 


BaH 2 


expt. 


XPS 


~5 a 


82 






XPS 


>2.2 6 


69 


LiH 


expt. 




4.99 


85 




calc. 


GWA 


5.24 


85 


LiD 


expt. 




5.04 


78 




calc. 


LDA 


2.84 


78 






GWA 


5.37 


78 



"Charging effects are taking into account in the experiments. 
''If there are no charging effects, the band gap would be twice as 
large. 



mirrors since in fully hydrogenated Y-Mg alloys Nagen- 
gast et al. found that fee YH 3 coexists with 7-MgH 2 i2& (3- 
MgH2 seems to be very different both structurally^ and 
optically. Its density is much larger and the calculated 
band gap for this material turns out to be considerably 
smaller than that of both a- and 7-MgH2 (see Table HV|l . 
The augmented plane wave (APW) calculations under- 
estimate the band gap considerably, the GGA value of 
2.35 eV (again by Vajeestonii) is much more reliable. 
The difference between the calculated band gap using 
GGA and the one measured for a-MgH2 is 1.4 eV. As- 
suming that the same scissors-operation can be applied 
we expect an experimental gap of 5.7 eV for 7-MgH2 and 
3.75 eV for /?-MgH 2 . 

We now compare MgH2 to related materials such as 
MgF 2 , MgO, MgS, MgSe, and other alkaline-earth and 
alkali hydrides. In Table E] the band gaps of these ma- 
terials are listed. For all of them it turns out that LDA 
underestimates the measured gap. Again, as for MgH2, 
the GW approximation seems to give a very good agree- 
ment between experiment and calculation for MgO, LiH 
and LiD. The same is true for the alkali halides J- MgH2 
and CaH2 seem to be very similar. Both materials are 
wide band gap insulators and have a valence band that 
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is predominantly determined by hydrogen orbitals^i^a 
For MgO and MgS both valence and conduction bands 
are determined by the anions (O or S)S as well and again 
the same holds for the alkali halides^ 

Since both the valence and conduction band of MgH2 
are formed by H states, the band gap of the alkaline- 
earth hydrides is expected to be almost independent of 
the metal as is the case for the alkali halidesp^I The 
remaining dependence results from the influence of the 
metal-ion on the lattice parameter, and the influence of 
hybridization of the conduction band states with metal- 
ion s- and d-states. From XPS experiments by Franzen et 
alml it seems that the onset of transitions in the valence 
band regions start for CaH2, Srfi2 and BaH2 all at about 
2.5 eV. Since charging effects are taking into account, this 
gives a gap of about 5 eV quite close to what we have 
found for a-MgH2. It is striking that LiH and LiD also 
have a gap of 5.0 eV. 

V. CONCLUSIONS 

In this study we use a novel experimental setup for op- 
tical transmission and ellipsometry measurements. This 
setup facilitates greatly measurements of the optical 
properties and dielectric function of metal hydrides in 
a hydrogen environment. It is possible to control the gas 
pressure from 1 mbar to 100 bar, and thus the composi- 
tion of metal hydrides over a wide range. The tempera- 
ture, pressure and resistivity are monitored in situ during 
hydrogenation of a sample. We determine the dielectric 
properties of a-MgH2 and PdH x and find that MgH2 is 
a transparent, colour neutral insulator with a band gap 



of 5.6 ± 0.1 eV. The transparency over the whole visible 
spectrum is ~80% (for a 100 nm thick film). The exper- 
imentally determined dielectric function in the photon 
energy range between 1 and 6.5 eV of a-MgH2 is in very 
good agreement with very recent calculations using the 
GW approximation. If we assume that calculations of 
the band gap for 7- and /3-MgH2 underestimate the ex- 
perimental gap by the same amount as in a-MgH2 we 
expect an experimental gap of 5.7 eV for 7-MgH2 and 
3.75 eV for /3-MgH 2 . 

In a coming publication we shall show that the dielec- 
tric function of MgH 2 determined here can be used to 
explain the large optical absorption (the so-called 'black 
state') of Mg-based alloys where metallic Mg and insu- 
lating MgH2 nanodomains coexist. 
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